Origin of the pseudo-gap and its influence on superconducting state 
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When holes move in the background of strong antiferromagnetic correlation, two spatial scales 
emerge which lead to a much reduced hopping integral with an additional phase factor. By taking 
these two effects into consideration, we propose an effective Hamiltonian to investigate pseudogap 
in cuprates. We argue that pseudogap is the consequence of dressed hole moving in the antiferro- 
magnetic background and has nothing to do with the superconductivity. In the normal state, The 
pseudogap opens near the antinodal region and vanishes around the nodal region. In the supercon- 
ducting state, the d-wave superconducting gap dominates the nodal region, while the pseudogap 
dominates the antinodal region. A two-gap scenario is concluded to describe the relation between 
the two gaps. 



PACS numbers: 74.72.-h, 74.72.Kf, 74.25.Jb, 74.20.-z 

I. INTRODUCTION 

One of the most fascinating properties of the cuprates 
is the opening of a gap (pseudogap) above T c in the 
underdoped and optimally doped regime, where most 
abnormalities observed^. The origin of the pseudo- 
gap and its relation to the superconducting (SC) gap 
are fundamental questions to realize the physics under- 
line the high-temperature superconductors. Although 
much progress has been made, the issue remains open. 
There are two distinct scenarios on the relation between 
the pseudogap and the SC gap£. 1) The one-gap sce- 
nario: The pseudogap is viewed as the precursor of 
the SC gap, reflecting pair fluctuation above T c , and 
would acquire phase coherence below T,A The argu- 
ment was based on angle-resolved photoemission spec- 
troscopy (ARPES)£A elect ron tunnelling 7 , and thermal 
transport measurements®, etc; 2) The two-gap scenario: 
The pseudogap is not directly related to the SC gap, but 
emerge from some ordered states such as antiferromag- 
netic (AF)i2, staggered flixsi 3 -, stripe^, spin or charge 
density wave 15 , and orbital circulating currents^, etc. 
and competes with SC gap. Models account for these or- 
dered states capture some aspects of the pseudogap, but a 
theoretical framework which gives a full picture from the 
pseudogap to SC state is not yet established. Neverthe- 
less, many experiments, including ARPES^, electronic 
Raman scattering^, and elastic neutron diffraction 11 , 
etc. seem to favor the two-gap scenario. More recently, 
the improved ARPES datai 7 - - — showed that the gap 
opens near the antinodal region, and vanishes near the 
nodal region in the pseudogap state, which leads to the 
arc structure of the Fermi surface. In the SC state, the 
pseudogap dominates the antinodal region, while the sim- 
ple d-wave type gap dominates the nodal region. 

In this Letter, an effective single particle Hamiltonian 



is proposed to study the pseudogap. The main idea came 
from the well known fact that moving holes in the anti- 
ferromagnetic background leads to two spatial scales: a 
long distance one and a short distance one. On the large 
spatial scale, the effective hopping integral is much re- 
duced due to accumulated strong AF correlation, while 
on the small spatial scale, the moving holes obtain an 
additional phase factor due to the surrounding AF corre- 
lation. Thus, pseudogap came directly from the AF cor- 
relation and is irrelevant to the superconductivity. Based 
on the effective Hamiltonian, we find that the pseudo- 
gap opens near the antinodal regime and vanishes around 
the nodal region. Its doping and momentum dependence 
qualitatively agree with recent ARPES measurements. 
The arc structure of the Fermi surface is also obtained, 
with its length expanding from underdoping to overdop- 
ing. In the SC state, the simple d-wave SC gap dom- 
inates the node region, while the pseudogap dominates 
the antinodal region. The resultant gap also agrees well 
with experiments. 

The neutron scattering^ and quantum oscillations 
measurements^ 2 , show that in the underdoped region, an- 
tiferromagnetic correlation remains on the Cu02 plane. 
Its existence has also been found numerically^ 3 .. For each 
spin, its four nearest neighbors compose a plaquette. As 
a hole is introduced into the center of this plaquette, it 
will sense an effective magnetic field originated from its 
four neighboring opposite spins. Correspondingly, there 
is a doping dependent gauge flux $ passing through the 
plaquette. In its journey of wanderings, the hole is af- 
fected by the staggered flux filed in which <f> and — $ ap- 
pear alternatively. Therefore, in the small spatial scale, 
when the hole hops, it encounters a phase shift 5/4 (or 
—(5/4) via Aharonov-Bohm effect, where S — 27r<&/ <&o 
($o = hc/e is the flux quanta). On the other hand, in 
the large spatial scale, the effective hopping, represented 
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TABLE I: The ED results of the spin-spin correlation function 
and hopping matrix in the 20-site cluster with different hole 
concentration. 

nh {Sj ■ Sj) h h ^3 

-0.3454 

1 -0.2745 0.0703 -0.0297 0.0125 

2 -0.2179 0.1483 -0.0350 0.0355 

3 -0.1670 0.2116 -0.0469 0.0510 

4 -0.1399 0.2635 -0.0594 0.0732 

5 -0.0855 0.3126 -0.0709 0.0960 



by the hopping matrix I n — (c^Cj + cfcj} n (n = 1, 2, 
and 3 for the nearest, second nearest, and third nearest 
neighbor respectively), evaluated by all possible configu- 
rations, is substantially reduced by the strong AF corre- 
lation. Take the t — J model as an example, we derive the 
following effective single particle Hamiltonian by consid- 
ering the effects of two spatial scales as, 

H = Y1 {^dt a e ka + he) +Y2 ek ( d L d k* + 4o e kv) ,(1) 

k<7 ka 

where ~f k = -2{h + J'xo)(e 4<5/4 cos k x + e~ iS / A cos k y ), 
and e k = — 4/ 2 cos k x cos k y — 2I 3 (cos 2k x + cos 2k y ) — [i, 
Xo = 2 r Zcr(^iCT^ff+ e io- e io-) i s the uniform bond order, J' = 
3J/8. The summation is restricted in the AF Brillouin 
zone. The effective Hamiltonian is diagonalized with the 
quasiparticle dispersion ef — e k ± |7fc|. 

To determine the phase factor (5/4, we calculate the 
effective staggered antiferromagnetic field, B e ff, which 
can be estimated by the spin-spin correlation between 
a spin 1/2 sits on the site of the given hole and its sur- 
rounding spins. That is B e ff = J J2jeNN(^i'^j)/9l J 'BS z 
with g, and \xb denoting Lande factor and Bohr mag- 
neton, respectively. Correspondingly the staggered flux 
$ = B e ffa 2 . a is the lattice constant, which is taken as 
0.45nm throughout this paper. {Si ■ Sj) is the nearest 
neighbor (NN) spin-spin correlation, which is calculated 
by the exact diagonalization calculations (ED) as shown 
in Table. HI Since the effect of long-range spin-spin corre- 
lations and spatial variation of the effective magnetic field 
are neglected, the resultant flux and corresponding phase 
factor are overestimated to some degrees. The effective 
hopping matrix /„ is also evaluated by the ED technique. 
For more hole cases, the hopping matrix should be di- 
vided by a factor approximately equal to hole number 
rift. In the above calculations, we have assumed that the 
extend t — J model can be used to describe the many 
body effect qualitatively, t (about OAeV in the real sys- 
tem) is taken as energy unit, t% — —0.3, i 2 = 0.2 and 
J = 0.4. Comparing with the slave-boson treatment* 2 *, 
the effective hopping terms are all enhanced due to many 
body effect. 

We further include the effect of the spin fluctuation 
under the random-phase approximation (RPA) 25 . Since 
partial spin fluctuation has been taken into account for 
the contribution of additional phase factor in Eq. (1), 



we introduce an adjustable parameter to set the AF in- 
stability at experimental value x ~ 0.03 to avoid the 
overestimating 2 —. The bare spin susceptibilities are ex- 
pressed as 

Xo(<7, w) = \ + F— ) + «-(*+- + 

k 

k 

with abbreviations F m > = /(£ *' ^f^l+J _ +,-). 

k-f-g k 

K± = 1 ± {Re^ k Re^ k+q + Im^ k Im^ k+q ) /\^ k+q \ and 
k 3 = (Rej k Imj k+q - Imj k Rej k+q )/\j k j k+q \. Since ± 
band is not independent due to the AF correlation, the 
nondiagonal term Xo(q,Q + Q,oj) arises from the umk- 
lapp processes. Under this treatment, the Dyson's equa- 
tion can be expressed as G{k,iLu n )~ l = Go(k,iuj n )~ l — 
T,(k,iu) n ) with £i(k,ioj n ) = Y. q , m J lx t { ( l,i UJ m)G(k - 
q,ioJ n — iojfn), where J q = J(cosq x + cosq y ), uj n and 
ijj m are the fermionic and bosonic Matsubara frequency. 
Both the Green's function and self-energy are 2x2 ma- 
trices. The spectral function is obtained via A(k,uj) = 
-(l/n)Im[G(k,u> + iS)]. 

The quasiparticle band shown in Fig. [1] is rather rigid 
against doping along the nodal line with much reduced 
bandwidth 0.6-0.7i(~240-280meV), well consistent with 
the ARPES measurement^. The upper and lower band 
coincide at (tt/2, 7r/2). This means that no full gap opens, 
unlike the Mott insulator state. On the other hand, the 
lower band around (ir, 0) shows clear flatness below Ef. 
Correspondingly, a pseudogap (denoted by Ape) opens 
around the antinodal regime. The distance between the 
flatness part and Fermi level decreases with increasing 
doping. The results obtained here qualitatively coincide 
with the ARPES data 26 . It should be noted that as far as 
the QP dispersion is concerned, even the effective single 
particle Hamiltonian (Eq. (1)) has already succeeded to 
describe the main features mentioned above as shown in 
the insert of Fig. [TJ 

Fig. [5] shows the evolution of the pseudogap in the mo- 
mentum space at different doping density. The magni- 
tude of pseudogap is determined by evaluating the mini- 
mal distance of the lower band from Fermi energy in the 
given direction. Apc{0) decreases with increasing 9, and 
disappears at a critical value of 9. This can be easily 
realized from Fig. [TJ where the gap opens near the antin- 
odal regime due to the flatness of QP dispersion induced 
by the staggered magnetic flux, while the gap closes due 
to the crossing of the QP dispersion over Ef. Such be- 
havior had been confirmed by various experiments^— . 
Both the pseudogap region and its magnitude decrease 
with doping as evidenced by ARPES data£. The largest 
Apg decreases from about 92meV at deep underdoped 
region to about AQmeV at optimal doping with almost 
linear doping dependence. The numerical results consis- 
tent with ARPES measurement, but the range of angle 
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FIG. 1: Lower band QP dispersion for different hole concen- 
tration under RPA treatments. Insert is the results obtained 
by direct diagonalizing the Eq. (I). 
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FIG. 2: Pseudogap as functions of angle 8 for different dop- 
ing. Solid lines, and open squares correspond to results with 
and without RPA treatment, respectively. Insert shows the 
experimental data on Bi2212 obtained from Ref. — , squares, 
circles, and triangles are for UD75K, UD92K, and OD86K, 
roughly corresponding to underdoped (x=0.1), optimal doped 
(x=0.15), and overdoped (x=0.2) case. The dotted line are 
guided for eyes. 



seem to have about 10% reduction^. Results obtained by 
direct diagonalizing the Eq. (1) show similar doping and 
momentum dependence. However, the magnitude and 
the variation range of angle are more reduced. Compar- 
ing with the experiments, the present pseudogap vanishes 
somewhat slower with increasing 9. In fact, the value of 
angle 9o(Apo{9) = 0) is intimately related to the AF 
correlation in momentum space, large 9q means smaller 
AF correlation. This implies that the present treatment 
is more appropriate to describe situations with strong AF 
correlation. In the present treatments, we do not adopt 
any pairing potential, which is distinct from the previ- 
ous slave-boson calculation^. Furthermore, the behavior 
of pseudogap differs from the early ARPES discovery^, 
and is far from that of the simple d-wave gap. Therefore, 
the pseudogap is not related to the SC pairing. Our the- 
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FIG. 3: Density plots for integration of spectral intensity 
times Fermi function over an energy interval [—0.025, 0.025]t 
around the Fermi energy for different doping density. 



oretical results imply that the pseudogap exists even in 
the overdoped range. 

The evolution of the Fermi surface (FS) is depicted in 
Fig. [3J The intensity is strong near the nodal region, 
and becomes weak gradually toward antinode due to the 
opening of the pseudogap in this region. Thus, an arc 
FS structure^! forms around the nodal region. Its length 
expands with increasing hole concentration. The arc ex- 
tends to a large hole FS above x — 0.20, where the pseu- 
dogap vanishes gradually. So, the arc structure is a direct 
consequence of the pseudogap. The evolution of the arc 
structure and its length also qualitatively agree with the 
recent ARPES experiment^. 

Now we discuss the existence and influence of pseu- 
dogap in the SC state. Based on experimental 
observations^— and above analysis, we start with 
the assumption that the total gap of the quasiparti- 
cle dispersion contains two components: the pseudo- 
gap, and the SC gap of the standard d-wave BCS form 
Afe = Asc{cosk x — cosky). After some manipulations, 
the corresponding Green's function is Gsc(k,iuj) = 
[G PG (k, iu)- 1 + AlGpc{k, —iuj)]^ 1 , where G PG {k,iu) 
contains pseudogap Apg(0), is the diagonal Green's 
function in the pseudogap state under RPA treatment. 
We then compare our results with the experiments to see 
if Apg(0)/Asc « 1 (one-gap scenario) or it is substan- 
tial (two-gap scenario). We take the x = 0.15 case, which 
is near the optimal doping, as an example. In the nodal 
region, the gap increases almost linearly with decreasing 
angle 9 as shown in Fig. [4j This is the typical d-wave be- 
havior, showing that d-wave SC pairing dominates this 
region. In fact, the gap comes from the SC condensation 
entirely because the pseudogap vanishes along the arc. 
In the antinodal region, large deviation from the stan- 
dard d-wave is rather obvious. The shape and the value 
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FIG. 4: The total gap in the SC state for different magnitude 
of Asc- The doping density is set to be x = 0.15, near the op- 
timal doping. Insert shows experimental data at optimal dop- 
ing, together with our theoretical results with Asc = 0.04. 
Circles are for LSCO (Ref. and triangles are for Bi2201 
(Ref. 

of the gap in the SC state with small value of Asc are 
quite similar to that of the pseudogap state. Even one 
increases Asc up to 0.06 such that 2A$c > ^pg(O), the 
shape keeps unchanged and the value of the gap enhances 
little. Therefore, the pseudogap dominates the antinodal 
region. Between the two regions, the pseudogap and the 
SC gap coexist and compete with each other, so the de- 
viation from the simple d-wave enhances gradually when 
approaching to the antinodal region. These features con- 
sist well with recent experiments on the optimal doped 
Bi2213±£, Bi220li2, and LSC02S, as shown in the inset 
of Fig. H 

The degree of the deviation depends on the ratio 
of Asc/Apc^O). Smaller ratio means larger devia- 
tion. Such effect has been manifested in optimal doped 
Bi2201 1? and Bi221222. The pseudogap in the two 
cuprates is almost the same at optimal doping, but the 
SC critical temperature in the former is three times of 
the latter. Correspondingly, the deviation from the sim- 
ple d-wave is much enhanced in Bi2201. The enhance- 



ment of the deviation with decreasing doping can also 
be explained following this way. Such deviation had also 
been obtained in the previous one-gap scenario by apply- 
ing the spin fluctuation theory^. However, the correc- 
tion was too small to account for the large deviation in 
underdoped Bi2201 and Bi2212. Based on our analysis, 
the two-gap scenario seems more appropriate to account 
for the observed gap behavior. The experiments sup- 
port the one-gap scenario, mostly came from the earlier 
measurements 4 . The data near the nodes were not easy 
to distinguish at that time. Additionally, some measure- 
ments, such as the thermal transport properties- may 
concern mainly on the nodal region, where the simple d- 
wave SC gap dominates. At last, we would like to point 
out that the enhancement of the SC gap with doping in 
the underdoped cuprates is an outcome of the decreasing 
pseudogap, which leads to the suppression of the density 
of state near the Fermi energy. 

In conclusion, we have proposed an effective single par- 
ticle Hamiltonian extracted from the exact diagonaliza- 
tion studies of the t — J model to investigate pseudogap 
in the high-Tc cuprates. The main idea is to take the ef- 
fects of two spatial scales, one leads to much reduced hop- 
ping and one adds to a phase factor, both reflect physics 
of doped holes moving in AF background, into account. 
Our results revealed that the pseudogap is the single par- 
ticle property, and is not directly related to the SC state. 
It opens near the antinodal region and vanishes around 
nodal region. The doping and momentum dependence 
of the pseudogap are qualitatively consistent with recent 
ARPES measurements. The arc structure of the Fermi 
surface is a natural product of the pseudogap. In the 
SC state, the simple d-wave SC gap dominates the nodal 
region, while the pseudogap dominates the antinodal re- 
gion. The two-gap scenario is therefore concluded. 
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